High-j single-particle neutron states outside the N=82 core 
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Abstract 

The behaviour of the ^13/2 and /19/2 single-neutron strength was studied with the (a,^He) reaction 
on "^^^Ba, ^^''Ce, ^''^Nd and ^^"^Sm targets at a beam energy of 51 MeV. The separation between the 
single-neutron states ^13/2 and /i9/2 was measured in A^=83 nuclei with changing proton number. 
To this end spectroscopic factors for states populated in high-^ transfer were extracted from the 
data. Some mixing of i=b and 6 strength was observed with states that are formed by coupling 
the /7/2 state to the 2^*" and 3^ vibrational states and the mixing matrix elements were found to 
be remarkably constant. The centroids of the strength indicate a systematic change in the energies 
of the ii3/2 and /i9/2 single- neutron states with increasing proton number that is in quantitative 
agreement with the effects expected from the tensor interaction. 
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The description of nuclear structure in terms of single nucleons moving within the nuclear 
mean field is at the root of much of the understanding of atomic nuclei. The properties of 
the available single-particle orbitals close to the Fermi surface are essential ingredients in the 
understanding of many phenomena, from the description of nuclear excitations within a shell 
model to the development of collectivity in low-lying vibrational states. There is increasing 
experimental evidence, emerging in situations where shell structure can be tracked across 
a range of neutron excess, that the sequence of single-particle states is not static and can 
change with neutron-to-proton ratio, indeed changing the 'closed shells' in some regions 
of lighter nuclei that are far from stability. For example, changes in single-particle states 
occurring in neutron-rich oxygen isotopes have been observed that destroy the familiar magic 
gap at = 20 close to the neutron dripline, whilst producing a new one at = 16 1|. 
In heavier nuclei, systematic shifts in energies of single-particle states have been observed, 
indicating similar evolution of shell structure with changing neutron excess j^. 

Recently, Otsuka and co-workers identified the tensor part of the nucleon-nucleon inter- 
action as a general feature that appears to drive these observed trends in the evolution of 
nuclear shell structure in different mass regions [3]. This is an interesting finding since, 
beyond the quadrupole moment of the deuteron, there have been only a few instances where 
this inherent component of the nucleon-nucleon force manifests itself in nuclear structure in 
an obvious fashion [3]. The tensor component of the nucleon-nucleon interaction has now 
been shown to have a systematic effect on single-particle energies, an effect that can even 
cause the making and breaking of magic numbers j^. Efforts are now underway to include 
tensor cornponents into Hartree-Fock (hf) mean-field models where they had been largely 
neglected S Q]- 

Care needs to be taken when comparing theoretical predictions of single-particle ener- 
gies with experimental data. Often little is known about the single-particle nature of the 
states in question. The energy of the lowest state of a given spin-parity can give misleading 
estimates for the energy of the corresponding single-particle orbital if there is significant 
fragmentation of strength. Nucleon transfer reactions are a good probe of such fragmen- 
tation. However, even if experimental spectroscopic factors have been determined using 
such methods, their absolute values are often of questionable meaning as they depend quite 
sensitively on choices made in the calculation of reaction cross sections. For example, small 
changes in the bound-state radius of the transferred particle, or differences in optical-model 
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potentials, can have a dramatic effect on the magnitude of the calculated cross sections 
without appreciably changing the predicted angular distributions. The comparison of iso- 
lated experiments on different targets can therefore be fraught with difficulty arising from 
the details of the reaction models employed, in addition to any systematic experimental 
differences. Careful studies across all targets of interest, using the same experimental meth- 
ods and the same choices in the extraction of spectroscopic factors, can alleviate most of 
these difficulties, at least in obtaining good relative spectroscopic factors. The kinematic 
characteristics of a particular transfer reaction also influence the reliability of spectroscopic 
factors. The assumptions inherent in the reaction models are best met where there is good 
momentum matching, ensuring relatively large cross sections for the direct transfer process. 
Many studies have been performed on reactions, such as (d,p), which are well matched for 
low-i transfer. However, there are fewer systematic studies where there is a large reaction Q 
value, and consequently higher momentum transfer, leading to large cross sections for direct 
high-£ transfers. Such reactions are therefore needed for reliable measurements of high-j 
orbit als. 

The work reported in Ref. 2f| describes the use of the {a,t) reaction on stable tin targets 
to populate high-j single-proton states in Sb isotopes. The relative spectroscopic factors for 
i=4 and 5 transfer were extracted using the same parameters in distorted-wave Born ap- 
proximation (dwba) calculations and a common normalisation factor. The results identifed 
the lowest 11/2^ and 7/2+ states with the dominant fragment of /iii/2 and (77/2 single-proton 
states, with strengths constant to within ±10%. No single potential could describe the en- 
ergetic separation of these states as a function of neutron number until the introduction of a 
tensor term in the residual interaction j^, ^ . The variation of the lowest known experimental 
states with the requisite spin and parity points to the possibility of a similar evolution of 
single-particle energies for neutrons outside N=82 in particular, for the 1^113/2 and i'hg/2 
orbitals. However, the single-particle nature of the states involved is somewhat uncertain as 
no previous study of had been undertaken using reactions that are well matched for high-£ 
transfers. 

Here we present a systematic study of spectroscopic factors across the stable metallic 
A^=82 isotopes using the (a,^He) reaction, which has good kinematic matching for £ ~ 6. 
We have identified states populated with high-£ transfer and extracted relative spectroscopic 
strengths for i=5 and 6, allowing the centroids of ^13/2 and hg/2 strength to be determined. 
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Since these states have no radial nodes, the comparison of their relative energies is largely 
insensitive to anything other than specific features of the residual interaction, as was the 
case for the states in Sb isotopes studied earlier. 

A beam of a particles with an energy of 51 MeV was supplied by the Yale tandem Van 
de Graaff accelerator and was used to bombard targets of ^^^Ba, ^^°Ce, ^^^Nd and ^^^Sm. 
Oxygen is an almost inevitable contaminant in such targets due to their chemical reactiv- 
ity. To avoid complex handling procedures associated with self-supporting metallic targets, 
material was evaporated from enriched oxides onto supporting carbon foils of thickness 20- 
40 /ig cm~^, which were then stored in an argon atmosphere with only brief exposure to air. 
Reactions on both oxygen and the carbon backing complicated the analysis, as discussed in 
detail below. The Yale Enge split-pole spectrometer was used to momentum analyse the 
products from the reactions. A gas-filled ionisation chamber and plastic scintillator, at the 
focal plane of the spectrometer, were used to isolate ^He ions from other reaction products 
on the basis of their energies and energy-loss characteristics. A Si monitor detector was 
placed in the target chamber at 30° for diagnostic purposes. The beam dose was measured 
using a Faraday cup and beam-current integrator. The beam currents were typically in the 
range 10-30 pnA. Data were collected using the same experimental set-up in two separate 
experimental runs. 

In order to extract absolute cross sections, the product of the target thickness and spec- 
trometer aperture was calibrated using elastic scattering at low beam energies. In the first 
run this was done with a spectrometer angle of 30° at an incident energy of 10 MeV. In the 
second run, elastic scattering of 20-MeV a particles was measured at 20°. The deviation 
from Rutherford scattering under both of these conditions is predicted to be less than 0.5% 
in calculations using several different optical-model parameters. The entrance aperture of 
the spectrometer was kept constant throughout each run to avoid any systematic differences. 
Using the nominal value of the entrance aperture, the ^^^Ba, ^^°Ce and ^^^Nd targets were 
found to be in the range 100-150 yUg cm~^; the ^^^Sm target was thinner at ~40 /ig cm~^. 
No significant variation in the normalisation from yield to cross section was observed be- 
tween measurements made at the beginning and end of each run, and between the two runs, 
showing that the targets are fairly robust. 

Typical momentum spectra for the observed ions are shown in Fig. [H where the main 
N=82 reaction strengths are via i=5 and 6 transfer, with the exception of some i=3 tran- 
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sitions leading to the 7/2~ ground state of the residual nucleus. An energy resolution of 
~70 keV was obtained. Both £=5 and £=6 strengths are split into two groups in each iso- 
tope. Also apparent in the spectra are intense lines arising from reactions on ^^O and ^^C 
target contaminants. These contaminant peaks obscure some of the lines of interest at cer- 
tain angles. However, the kinematic characteristics are such that the momenta of the ions 
from contaminants change more quickly with angle than for the reactions on N=82 targets, 
allowing angles to be found where each peak of interest is free from contamination. Weaker 
contaminant transitions from ^^Si(a,^He) reactions were also present in some spectra, par- 
ticularly those from the ^^"^Sm target. They were easily identified using known excitation 
energies and from the momentum variation with angle and were not close to any peaks of 
interest . 

Data were taken at laboratory angles of 6°, 11° and 20° for all targets. The first angle 
was chosen to be as close to the peak yields for i=5 and 6 as possible without resulting in 
count-rate problems, whilst 11° and 20° are angles sensitive to the differences in the shapes 
of angular distributions. For the ^^^Ba target, where contaminant ion groups cause most 
problems, spectra were also taken at 30°. Since spectroscopic factors are most reliable when 
measured at the peak of the angular distribution at forward angles, measurements were also 
made at 30° using the ^^''Ce target as a consistency check on the reliability of spectroscopic 
factors obtained at the larger angles. 

Examples of angular distributions are shown in Fig. [2], together with dwba calculations 
using optical-model and bound-state parameters listed in Table I. The optical potentials were 
successful in reproducing angular distributions of high-£ transitions in Sn(a,t) reactions M - 
The exact finite-range calculations were performed using the computer code ptolemy [12| . 
The comparison with data indicates the quality of the dwba calculations in reproducing 
;he angular variation in cross section. Spin assignments were taken from the literature 
although, in general, the measured angular distributions agree with these 
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spin assignments on the assumption that i=5 and 6 transitions lead to the population of 
9/2^ and 13/2+ states, respectively. No strong transitions were observed with transfers 
other than i=3, 5 and 6. 

In ^^^Ce there is a doublet composed of a 13/2+ state at 1354.52(9) keV and a 9/2~ 
state at 1368.7(2) keV 1^] that are unresolved in the momentum spectra. Yields for these 
two states were extracted by fitting two Gaussian functions to the data with positions fixed 
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using the known excitation energies. The widths were fixed to values reproducing the shape 
of neighbouring singlet peaks, such that the peak heights were the only free parameter in 
the fit. This procedure was successful at 20° and 30°. At 6° and 11° satisfactory fits were 
not obtained due to an increasing background in the spectra from inelastically scattered a 
particles that cannot be completely separated from ejectile '^He ions at forward angles. An 
alternative method, fitting two calculated £=5 and 6 angular distributions to the summed 
yield of the doublet led to ambiguous results, with a long fiat trough in the space, due 
to the similarity in shape of the two distributions. 

Table II gives the differential cross sections measured in the current work, where ab- 
solute numbers typically have an error of ~ 7% whilst relative numbers are good to 5%. 
Where a cross section value is not listed in the table, that particular angle was affected 
by contamination by oxygen or carbon ion groups. In addition, data were not taken at 
30° for ^^^Nd and ^^^Sm targets and yields could not be extracted from the ^^^Ce doublet 
at forward angles, as discussed above. Also shown are the spectroscopic factors obtained 
using a single normalisation for both £ = 5 and 6 transfers across all isotopes and transi- 
tions listed. The spectroscopic factors at different angles typically varied by less than 10% 
from the average values listed in Table II, apart from those arising from the Ce doublet 
peak where the observed variation was less than the larger errors arising from the fitting 
procedure. The relative values of the spectroscopic factors are subject to an uncertainty of 
~ 10%, based on the error in the relative cross sections and analyses performed using a va- 
riety of different sets of optical-model parameters. Spectroscopic factors were not extracted 
for the I = ?) transitions to the residual ground states; dwba calculations were not able to 
reliably predict the shape of these poorly-matched transfers. The spectroscopic factors show 
significant deviations from those extracted previously using (d,p) reactions 
and references therein), not only between different targets, but even for a given target with 
the same £-value, as might be expected when comparing spectroscopic strengths extracted 
from poorly-matched and well-matched transitions. The totals for I =5 and for 6 strengths, 
shown in Table III, are constant to within ±13% across the isotopes studied. 

High-£ transfers populate two states for both £ = 5 and 6 in each final nucleus; no 
other transitions were observed with these i transfers to a limit of 5% of the strongest 
transition. This fragmentation is likely to be due to coupling to vibrational excitations of 
the core, leading to predominately two-state mixing between O+j.^, ® 2^^13/2 and "i^axc ® 
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configurations in the case of tlie positive-parity states, and 0^j,g ® z^/i9/2 and (g) Z//7/2 
configurations for tlie negative-parity states. Simple two-level mixing calculations have been 
used to extract the mixing matrix elements from the measured spectroscopic factors and 
excitation energies. The values of the matrix elements are shown in Table III. They are 
remarkably constant across all the nuclei studied, leaving the extent of the mixing to be 
dictated by the energetic proximity of the unperturbed states, illustrated in Fig. [3] (a) 
and (b). More detailed calculations with particle-core coupling models have been done in 
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lere parameters were fitted to the available energies and spectroscopic factors 



20| . The results of these calculations are broadly in line with the experimental 



results of the current work, although they have been fitted to spectroscopic factors from 
(d,p) reactions so a detailed comparison is not appropriate here. 

The trend in the ^13/2 and hg/2 centroids is similar to that of the lowest 9/2~ and 13/2"'" 
states, as illustrated by the energy differences shown in Fig[3](c). However, the fall in energy 
with increasing Z is not as large and the crossing of the two lowest states is not observed in 
the single-particle energies. 

In the Z=51 isotopes, the observed evolution of the 7r/iii/2 — 7^97/2 gap as a function of 
neutron number, which is somewhat similar to the trends uncovered here, has been success- 
fully reproduced by the introduction of the tensor-force terms in shell-model and mean-field 
calculations j^, ^. Qualitatively, these shifts have a fairly simple effect: nucleons filling a 
j = i ± s orbit have an attractive effect with nucleons in j = i ^ s orbitals, whereas a 
repulsive contribution to binding energy arises from the interaction between nucleons where 
both are either in j=i + s ot j = i — s states . 

In the case of neutron states outside the A^=82 core, the tensor force will also produce such 
a shift arising from changes in the occupancy of various proton orbits. Experimental proton 
occupancies in N=82 nuclei are available from previous systematic measurements of proton- 
transfer reactions [2l|; over the stable isotopes, the nlgj/2 and vr2(i5/2 orbits were observed to 
fill at approximately the same rate. These proton orbitals have opposite effects on neutron 
binding energies, having different spin-orbit couplings, j = i — s and j = i + s. But the effect 
of the 7rl(77/2 state should dominate, due to the considerably larger spatial overlap with the 
nodeless radial wavefunctions of the high-j neutron states. An increased binding energy 
for zi3/2 neutrons, and a decrease for hQ/2, should therefore result from a tensor interaction 
with (77/2 protons. As the occupancy of this proton orbital increases, the two neutron states 
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should drift towards each other in energy. This quahtative expectation is in agreement with 
the current experimental results, as illustrated in Fig. [31 More quantitatively, the effect of a 
tensor force based on vr and p meson potential (discussed in detail in Ref. fs}]) is to reduce the 
difference between and /i9/2 neutrons by 0.180 MeV per additional proton in the 
orbital, using calculations of the relevant matrix elements using yl=140 harmonic oscillator 
wave functions [2^- In line with the qualitative expectations discussed above, the interaction 
with a proton is calculated to increase the gap by only 0.04 MeV. Combining these 

tensor matrix elements with measured proton occupancies 2l|, the relative changes produced 
by the tensor interaction on the ^11^/2 — 1^^9/2 gap can be deduced. A comparison of this 
calculation with the experimental single-particle energy differences is shown in Fig. [21(d), 
displaying reasonable agreement with the measured changes to the i'iiz/2 ~ ^^^9/2 gap. 

Calculations with mean-field models are more difficult in this case, compared to the 
situation with proton orbitals outside Z=50, due to the need to reproduce the details of 
the filling of protons in the close- lying TTg'i/2 and 7[d^/2 orbitals near the Fermi surface. 
These states are known experimentally to lie only ~100-300 keV apart 2l|. Recently HF 
plus BCS calculations have been published for neutron states outside N=82 [3], based on 
a modification to the SLy5 Skyrme interaction that introduces elements from the tensor 
force. It predicts a different pattern of filling proton orbitals compared to that deduced 
experimentally, namely sequential filling of the 77(77/2 level, followed by the 77^5/2 orbital. As 
a result, there are discrepancies between the calculated changes to the 1^113/2 — z^/?-9/2 gap as 
a function of Z and the experimental data presented here. 

No quantitative information exists on the single-neutron nature of high-j states in A^=83 
nuclei, other than those accessible in neutron-adding reactions on the stable metallic species. 
There are no 13/2+ states known in ^^'^ Xe; measurements in normal kinematics for a ^^^Xe 
gas target are compromised by reactions on gas-cell windows leading to contamination prob- 
lems. In lighter A^=83 nuclides, direct experimental information about low- lying 13/2"*" 
states is missing. However, the energies of both the lowest 13/2+ and 9/2^ states are known 
for Z > 62, even though spectroscopic factors are unavailable. Beyond Z = 62 there is a sud- 
den change in the energy systematics; the two states start to diverge in energy quite strongly 
(see Fig.[3]^c)). In these higher Z systems, the 7r/iii/2 state should begin to fill, with dramatic 
consequences for high-j neutron states due to their large radial overlap and an increasing 
repulsive contribution from a tensor interaction, amounting to 0.16 MeV per additional 
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proton [22] . However, any strong conclusions should be treated with some caution since 
the single-particle content of the lowest-lying states will not be known until nucleon-adding 
reactions are performed on radioactive beams. Indeed, it is towards the Z=64 sub-shell 
closure that the 3^ core excitation is expected to fall closest in energy to the centroid of 
^13/2 strength (see Fig. EJ^b)). Assuming the mixing matrix elements are of similar strength 
to those extracted here, this will lead to large mixing and significant fragmentation of i=6 
strength between two 13/2+ states. The energy of the lowest 13/2"^ state will therefore be 
a particularly unreliable estimate of the 213/2 single-neutron energy. 

In conclusion, spectroscopic factors have been measured for high-j states outside stable 
A^=82 isotopes revealing a systematic drift in ii3/2 and hQ/2 single-neutron energies toward 
each other with increasing proton number. There is quantitative agreement between the 
experimental trends and the shifts in single-particle energies predicted with a tensor inter- 
action. 

We are indebted to Takaharu Otsuka for providing us with the relevant tensor matrix 
elements and for many very fruitful discussions. We would like to acknowledge meticu- 
lous target manufacture by John Greene (Argonne National Laboratory) and Paul Morrall 
(Daresbury Laboratory). Peter Parker has provided invaluable assistance and advice, with- 
out which these measurements would not have been possible. We would also like to thank 
the staff at Yale University, in particular, their efficient facility operation. This work was 
supported by the UK Science and Technology Facilities Council and the US Department 
of Energy, Office of Nuclear Physics, under Contract Numbers DE-FG02-91ER-40609 and 
DE-AC02-06CH11357. 
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FIG. 1: Spectra from the (a,^He) reaction as a function of excitation energy taken at 20° in 
the laboratory frame, labehed by residual nucleus. Peaks are labelled with the i transfer of the 
transition. The large peak moving up through the spectra with increasing atomic number is due 
to reactions carbon and oxygen contaminants in the target. Small, broader peaks, particularly 
evident in the ^^^Sm spectrum, arise from silicon contamination. 

FIG. 2: Differential cross sections from the ^'^^Sm(Q,^He)^^^Sm reaction shown for transitions 
leading to the following states in the residual system: (a) 13/2+ state at 1105 keV, (b) 9/2~ state 
at 1423 keV, (c) 13/2^^ state at 2670 keV and (d) 9/2~ state at 1780 keV. Square points indicate 
i = 5 transitions, circular points i = 6 assignments. The lines are dwba calculations normalised 
to the data, for i=5 (continuous line) and for i=6 (dashed line). 

TABLE I: Optical-model and bound-state potentials used in the dwba calculations. Depths 
of the bound-state potentials, marked with were varied to reproduce the appropriate binding 
energies. 
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FIG. 3: (a) Experimental centroid energies of the neutron hg/2 single-particle strength for = 83 
nuclei (solid squares) compared to the quadrupole vibrational excitation energies of the correspond- 
ing even-even core (open symbols), (b) A similar plot for the ^13/2 centroid and octupole vibrations 
of the core. The difference in experimental centroid energies of 1^^13/2 and i'hg/2 strength (solid 
triangles) compared to: (c) the energy differences between the lowest 13/2+ and 9/2^ states (open 
triangles), and (d) predictions of a vr -|- p tensor interaction (grey band). The width of the grey 
band is generated by experimental uncertainty in the proton occupancies [21]. 



11 



TABLE II: Differential cross sections and spectroscopic factors of states populated by large i 
transfer in the (a,^He) reaction on stable N=82 targets. Excitation energies are taken from 
Ref. 3, 14, 15, 161. 
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TABLE III: Summed single-neutron strengths, energy centroids and mixing matrix elements with 
the state of corresponding formed by coupling the 1^/7/2 state to the 2+ or 3^ core vibrations. 
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